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Abstract 

During the runs of the PS 179 experiment at LEAR of CERN, we 
photographed an event of antiproton-iVe absorption, with a complete 
7r + — >fJ- + — >e + chain. From the vertex of the reaction a very slow en- 
ergy 7r + was emitted. The ir + decays into a [i + and subsequently the 
/i + decays into a positron. At the first decay vertex a muon neutrino 
was emitted and at the second decay vertex an electron neutrino and 
a muon antineutrino. Measuring the pion and muon tracks and ap- 
plying the momentum and energy conservation and using a classical 
statistical interval estimator, we obtained an experimental upper limit 
for the muon neutrino mass: m u < 2.2 MeV at a 90% confidence level. 
A statistical analysis has been performed of the factors contributing 
to the square value of the neutrino mass limit. 

PACS. 13.20 Cz, 14.60 Pq, 13.25 Cq 
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1 Introduction 

Determination of the absolute values of neutrino masses represents a most 
difficult problem from an experimental point of view. Evidence in favour of 
non-zero neutrino masses and oscillations obtained in most of the relevant 
neutrino experiments has made the physics of massive neutrinos a frontier 
field of research in particle physics and astrophysics. All the existing terres- 
trial and astrophysical data indicate that the neutrino masses are by many 
orders of magnitude smaller than those of other experimentally measured 
lepton and hadron masses. Such a low value is the most relevant reason 
for it being extremely difficult to extract the values of neutrino masses from 
experimental measurements, either directly or indirectly. This is the main 
reason why most experimental papers only report a confidence limit inter- 
val for such values. In this paper we report an upper limit for the muon 
neutrino mass obtained by measuring the radii of curvature of the pion and 
muon tracks of an event, in which a pion decays into a muon and a muon 
neutrino, recorded during the runs of the PS 179 experiment at the beam of 
antiprotons of LEAR at CERN. 

2 Experimental apparatus 

The PS 179 experimental apparatus was designed and built for the study 
of antiproton interactions with light and medium-light nuclei at the LEAR 
facility of CERN. The aim of the research was an experimental study of the 
interaction of antiprotons with nuclei at low energies. The results obtained 
provided information on the fundamental nucleon-antinucleon forces, on the 
interaction of antiprotons with clusters of bound nucleons, on the distribu- 
tion of nuclear matter, and on properties of highly excited nuclear matter, as 
well as a restriction on the possible amount of antimatter present in the early 
Universe [1]. In order to make the most of the information available on all 
the secondary charged particles produced in the reactions, in designing the 
experiment a choice was made of the visualization detection technique, and 
a self-shunted streamer chamber placed in a magnetic field was used. Such 
a detector had many advantages. It was a low density gas target, and, at 
the same time, it was triggerable, offering a 4n acceptance in which highly 
luminous localized particle tracks could be obtained [2] . The low density of 
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Figure 1: Layout of the PS 179 experimental apparatus, including: EM- 
electromagnet; SC-streamer chamber; HVPG-high voltage generator; TP- 
travelling platform; ES-electrostatic screening; VFCi_ 2 -wire chambers; Ci_ 6 
-scintillation counters; M!_ 5 -thin walls. 



the target medium, i.e. of the streamer chamber filling gas at 1 atm pres- 
sure, allowed to reveal long charged particle ranges and nuclear fragments. 
The experimental apparatus used is sketched in Figure 1. The detector pro- 
vided stereo pictures of the sensitive chamber volume. The details of the 
experimental apparatus have been given in ref. [3] where a more complete 
description of the characteristics and the performances of all the components 
of the setup can be found. 



3 Analysis of the measurements 

The runs of the PS 179 experiment were carried out with the streamer cham- 
ber filled at 1 atm with different gaseous targets: 3 He, 4 He, Ne. The stereo 
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pictures of events, recorded on films during the runs, were reprojected onto 
measuring tables for visual scanning. The events were measured with a dig- 
itized coordinatometer directly on the scanning tables. 

In one of the photographs of the iVe run, the beautiful event, reproduced 
in Figure 2, was observed. It represents a multi-nucleon annihilation [4] of 
an antiproton with a Ne nucleus (see, also, Figure 11 of ref. [3]). From 
the vertex of the interaction there "evaporated" a (1.98±0.02) MeV positive 
pion. From the two vertices of the pion and muon decays, three neutrinos 
were emitted: from the first a muon neutrino and from the second an electron 
neutrino and a muon antineutrino. The peculiarity of the event consists in 
the noticeable length of the 7r + and tracks and in their coplanarity. 

The temperature of the streamer chamber filling gas was 289 K, the pres- 
sure atmospheric, and the target density (0.80±0.01) mg/cm 3 . The magnetic 
field was 0.8 T (AB/B = 1(T 5 ). 

The tracks have been measured at JINR with a microscope- digitizer. The 
recent acquisition of such a microscope made it possible to newly measure 
the TTfie event, reconstructing and digitizing tracks with a very high accuracy. 
The pion track has been digitized at 347 points, the muon track at 2037 
points and the positron track at 180 points. All the measured points are 
approximately equidistant. Each point has been measured with a precision 
of about 10~ 2 mm. The coordinates of all track points have been used to 
deduce the radii of curvature and the ranges. Each particle transfers energy 
by ionization to the surrounding medium and slows down along its path. A 
Fortran code has been written to estimate the radii of curvature varying along 
the tracks and the total track lengths. For each track the best-fit circles were 
calculated taking into account a fixed number N of points. Starting from the 
first measured point and taking the subsequent N-l points, the first radius 
was calculated. If n is the total number of points measured along a track, 
shifting the N points by one along the track, the subsequent n-(N-l) radii 
have been calculated. We have chosen N = 25 (a value statistically significant 
in order to ensure fit stability) for each track. Program CIRCLE of the CERN 
library [5] has been included into the code to calculate the best-fit radii of 
curvature along the tracks. 

The initial pion momentum was estimated from all measured points along 
the first 10 cm of its track, since the range-energy relations (and tables) for 
neon at NTP reveal the energy loss of a 2.00 MeV pion to be negligible 
for path lengths not exceeding 10 cm. The average radius of curvature of 
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the first 10 cm of the pion track, obtained best-fitting the track with a cir- 
cle, is r = (9.79±0.09) cm corresponding to a momentum (p = 300Br) of 
(23.50±0.15) MeV/c. The agreement between the value of the initial pion 
momentum and its measured path range of (471.5±0.8) mm, shows that the 
pion decay occurred with a momentum of (0.05±0.10) MeV/c. The same 
parameters have been measured for the muon track. The initial muon mo- 
mentum was estimated from all measured points along the first 12 cm of 
its track, since, also in this case, the range-energy relations (and tables) for 
neon at NTP reveal the energy loss of a 5.00 MeV muon to be negligible for 
path lengths not exceeding 12 cm. The value of the average radius of curva- 
ture of the first 12 cm of the muon track, obtained with a best-fit circle, is 
r = (12.50±0.08) cm, corresponding to a momentum of (29.90±0.15) MeV/c. 
In this case, also, the agreement between the initial muon momentum and 
its measured path range of (2616. 0±0. 8) mm shows that the muon decayed 
with a momentum of (0.6±0.6) MeV/c. The positron was emitted with a 
momentum of (46.80±0.08) MeV/c. The angle between the tangent lines of 
the pion and muon trajectories at the decay vertex was (163. Oil. 0) deg, and 
the angle between the muon and positron tracks was (107. Oil. 0) deg. It 
should be pointed out that the njj,e event of Figure 2 is the only existing of 
this kind whose parameters could be all measured with high precision. 

3.1 Kinematic and dynamic constraints 

For kinematic and dynamic analysis of the reaction, the relations of conser- 
vation of momentum and energy have been used. The pion decay occurs in a 
plane. At the decay vertex, the muon and the neutrino have flight directions 
at angles <p and 9 to the pion flight direction, respectively. For momentum 
conservation 



Ptt — P/iCoscj) + p v cos0 
= p^sincf) i p u sin9. 



(1) 
(2) 




(3) 



vl = Pi + p\ - 2p 7r p fl coscf) 



(4) 
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Figure 2: PS 179 picture of "The Three Neutrino Event". An antiproton 
annihilates with the Ne nucleus. From the interaction vertex a n + of energy 
inferior to 2.0 MeV is emitted. The 7r + decays into a and subsequently 
the jd + decays into a positron. At the first decay vertex, a muon neutrino 
is, also, emitted, and at the second decay vertex an electron neutrino and a 
muon antineutrino. 
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and 



E n - E ll = \JpI+pI~ 2p~p^cos4> + ml 



(5) 



from which one can obtain the square of the neutrino mass 



m 2 = (E w - E^f - {pi +pl~ 2p 7r p fM cos<p) 
= (E n Ejj_ Pv )(E n -E IX + p v ) = 
= m*m** 



(6) 



where only the m* factor can approach zero and where m** » m* and 
always positive. It is clear that from relation (6), given the masses of the 
pion and muon, the momenta of the pion and muon, and the angle (p, one 
can calculate the square of the neutrino mass as a function of the above- 
mentioned parameters. 

3.2 On the upper limit of m v 

A point or an interval estimator of an unknown parameter is any statistic 
whose value is a meaningful guess for the value of the unknown parameter 
which is assumed to have some fixed value [6]. Sometimes experimental 
measurement may yield non-physical values, when the parameter value is 
near zero [7] . In these statistical procedure can be used for estimation 

of the upper limit value of the parameter. In the present analysis the neutrino 
mass is the parameter to be estimated. The upper limit of the mass of 
the muon neutrino [8] has a value near zero, and is about seven orders of 
magnitude smaller than the values of the measured experimental quantities 
of the 7T— ^[iv^ decay. In determining such a mass in an unbiased approach 
one would necessarily expect a non-physical result for m 2 to occur half of the 
time. This happens because of the uncertainties of the experimental values 
of the physical quantities. One should then make statistical inferences from 
observation of the squared mass m 2 ± Am 2 when m 2 is negative or near 
the non-physical region. Am 2 is the standard deviation. It is possible to 
choose a "classical" confidence level p for the squared mass such that the 
corresponding classical confidence limit m 2 cl is in the physical region. When 
the variable m 2 has a Gaussian distribution the upper limit is 



™>l,d = r m 2 + ZpAm 2 



(7) 
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at a p% confidence level, where = 1.036, Zqq = 1.282, Z95 = 1.645 
and Z97.5 = 1.960 (see ref. [9]). The classical confidence limit satisfies the 
probability statement 

P(m 2 Q < ml d ) = p. (8) 

This statement says that m^ ci has the probability p of being larger than the 
true value (m^), whatever m 2 really is. It means that, if one repeats the 
experiment many times, and each time one recalculates the value of m^ cl , 
then in p% of the cases m 2 vd will be greater than the true value m\. 

3.3 Monte Carlo events 

To obtain an interval estimate for the muon neutrino mass from our mea- 
surement of the scattering angle and of the momenta of the two charged 
particles of the decay, a statistical procedure has been applied. It takes into 
account the measured pion and muon momenta and the angle between them. 
These are distributed with Gaussian p.d.f.'s with means equal to the mea- 
sured values and standard deviations equal to the corresponding measured 
uncertainties. A code has been written. Using the Monte Carlo method 
a set of 10 5 decays has been generated. Each element of the set has been 
obtained extracting at random the momenta of the pion and of the muon, 
and the angle between the two particles assuming the mass values of m n = 
(139.57018±0. 00035) MeV and m M = (105.658369±0. 000009) MeV, as given 
in ref. [8]. The three variables were extracted independently of each other. 
Other relevant quantities were then calculated. The distributions of ml, of 
m* and of m** in relation (6) have been obtained. 

3.4 Analysis of the results 

Figure 3 shows the distributions of ml (upper part), m* (central part) and 
of m** (lower part). The ml distribution is Gaussian, with a negative mean 
value and relatively large standard deviation that depends on the uncer- 
tainties in the experimentally measured quantities. The m* distribution is 
Gaussian with a negative mean value and relatively small standard deviation. 
Both distributions have ratios of the standard deviation to the mean value 
of the same order of magnitude. The m** distribution is Gaussian but in the 
physical region, with a mean value of (59.7±0.2) MeV, about two orders of 
magnitude greater than m*. 
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Figure 3: Distribution of ml (a), of m* (b), and of m** (c) generated with the 
Monte Carlo method. Note the different scales along the horizontal axes. The 
data have been best fitted (dashed line) with Gaussians with the following 
respective mean values and standard deviations: (—11.1; 12.5) MeV 2 (a); 
(-0.186; 0.211) MeV (b); (59.7; 0.2) MeV (c). 
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The ml distribution leads to the squared muon neutrino mass 

ml = (-11.1 ± 12.5) MeV 2 (9) 

which is compatible with zero. According to the classical method [7], de- 
scribed in Section 3.2, equation (9) corresponds to the muon neutrino mass 
upper limit 

m v < V-ll.l + 1-282 • 12.5 MeV = 2.2 MeV (10) 

at the 90% confidence level. This is the lowest estimate of the upper limit 
value of the muon neutrino mass obtained with a visualizing detector, using 
all the directly measured kinematic and dynamic parameters of the n— >\xvy, 
decay, that is p w , p^ and <fi (see ref. [10] and [11] and references quoted 
therein). The last measurement of this kind of events has been performed 
with a helium filled bubble chamber [11]. Because of the density of liquid 
helium, the muon tracks were only 1 cm long, and composed by about (15±2) 
bubbles. The authors assumed that all the stopped pions decayed at rest. 
For this reason they did not use in the analysis the values of the angle <fi and 
of the momentum of each pion. It is quite likely that most of the low energy 
pions decayed in flight. 

Since then many experiments have been devoted to measuring the muon 
momentum with higher accuracy. Recently, Assamagan et al. [12], using a 
magnetic spectrometer equipped with a silicon microstrip detector, studied 
the decay of pions stopped in a graphite target. The authors measured for 
the muons an average momentum of (29. 79200±0. 00011) MeV/c. For the 
pions immediately before their decay, the estimated average kinetic energy 
was T^ = (0.425±0.016) eV, so most of the pions actually decayed in flight. 
Setting the energy of the pions and the angle <f> to zero, the authors deduced a 
squared muon neutrino mass of (— 0.016±0.023) MeV 2 and, according to the 
Bayesian approach, deduced the corresponding neutrino mass upper limit of 
m, < 0.17 MeV (C.L. = 0.9). 

Applying the procedure described in paragraph 3.2, we have estimated 
the upper limit of the muon neutrino mass from the data of ref. [12], taking 
into account the measured average energy value of the decayed pions and 
assuming a non-zero angular uncertainty A0 of 0.2 deg (differently from 
ref. [12]). In this case the distribution of m\ is totally in the physical region. 
The value obtained for the muon neutrino mass is m v = (800±80) keV. 
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If we set the momentum of the decayed pion equal to zero, as in ref. [12], if 
we assume the uncertainty of this value to be of the same order of magnitude 
as that measured for the muon, and take A0 = 0.2 deg, the distribution of ml 
extends in the non-physical region. The calculated upper limit of the muon 
neutrino mass is m u < 110 keV (C.L. = 0.9) using the classical approach [7]. 

These exercises reveal the necessity, for an accurate estimation of the 
squared muon neutrino mass, of precise and direct measurements of all kine- 
matic parameters of individual 7r— ^/xi/^ decay events, like those performed 
for the event reproduced in Figure 2. 

Analysing the distributions shown in Figure 3, we studied the character- 
istics of the statistical variables ml, m* and m**. Although both m* and m** 
have the same Gaussian distribution (being composed of the same physical 
quantities), the m** distribution is shifted toward large positive values. The 
m* distribution is close to zero and the m* values may be negative (this is 
true, also, for the ml distribution). 

Using the formulae of paragraph 3.1 we calculated the kinematic param- 
eters of the decay n— ^/x^, assuming the mass values of and m M as given 
in paragraph 3.3, and a value for the muon neutrino mass of 0.300 keV. The 
choice of this value for the neutrino mass is related to the actual precision of 
the pion and muon masses [8]. The mean value of the pion momentum was 
set to 0.05 MeV/c and of the angle (p to 163.0 deg in order to allow a di- 
rect comparison with our measurements. Applying the procedure described 
in paragraph 3.3 and running the Monte Carlo code with different values 
of the uncertainties in the input variables (within 10" 10 -hl0 MeV/c for the 
momenta p n and p^, and 10~ 5 -=-l deg for the angle 0), sets of 10 5 n— >\xv^ 
decays have been generated. Each decay has been obtained extracting at 
random, independently, the three Gaussian variables p n , p^ and <fi. The pro- 
cedure gives in output all the relevant kinematic and dynamic quantities of 
the decay, in particular, m* , m** and m^. 

In Table 1 we present the results for different sets of Monte Carlo events, 
all obtained with an uncertainty of 1° in the measurement of the (p angle. 
The choice of this value was again dictated by the uncertainty of our mea- 
surement. In the first two columns the uncertainties of the input variables 
are reported. The third, the fourth and the fifth columns show the mean val- 
ues and the standard deviations of the m*, m** and m 2 u distributions. Both 
m* and m 2 v distributions have negative mean values. The m** mean values 
are always positive. So, for all the experimental uncertainties reported in 
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Table 1: Uncertainties, mean values, standard deviations of the m*, m** 
and ml distributions. The data have been calculated assuming m u = 0.300 
keV. For momentum uncertainties less than 10~ 6 MeV/c the values remain 
constant, as shown in Figure 4. 

Table 1, the distributions of m* and w? v extend into the non-physical region. 
Using the same classical method applied to deduce the muon neutrino upper 
limit mass from the m 2 v distribution [7], we calculated an upper limit from 
the m* distribution. Figure 4 shows the muon neutrino mass upper limits 
calculated from the m 2 v distributions as a function of the uncertainties of the 
momenta of the 7r and \x for different values of the <fi angle, and for m v = 0.300 
keV. In the same Figure 4, the m* and m** (C.L. = 0.9) upper limits are also 
reported. The values calculated from formula (7) follow straight lines. The 
lines (a) and (b) have been drawn to guide the eye, the line (c) represents 
a constant upper limit value of 60 MeV. When the momentum uncertainties 
are larger than about 10 MeV/c, the lines (a) and (b) intersect at the upper 
limit of m v ~ 60 MeV. For momentum uncertainties less than 10 -9 MeV/c 
and for angular uncertainties less than 10~ 5 deg, the upper limit values of 
m v reach the assumed value of 0.300 keV, confirming the validity of the clas- 
sical approach [7] we used. Moreover, one can see that the limit values of 
m u attainable with present modern techniques correspond to measurement 
uncertainties in momenta higher than 10~ 4 MeV/c. 

To analyse the case of pions decaying in flight with higher momenta, we 
calculated the kinematic parameters of the ir^fiv^ decay for a pion mo men- 
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Figure 4: Behaviour of the upper limits of m v (open triangles, squares, stars, 
circles), m* (full triangles, squares, stars, circles), and m** (open rhombs) 
calculated from the ml, m* and m** distributions at 90% confidence level, 
as a function of the experimental uncertainties in the momenta, and for 
uncertainties in the angle <fi of 1.0 deg (open triangles, full triangles), of 10~ 2 
deg (open squares, full squares), of 10 -4 deg (open stars, full stars) and of 
10~ 5 deg (open circles, full circles). The lines (a) and (b) are drawn to guide 
the eye for A</>=10~ 5 deg and a p >10~ 9 MeV/c. The line (c) represents the 
upper limit value of 60 MeV. The values of 2.2 MeV and 85 keV are from 
Table 2. 
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turn of 200 MeV/c and neutrino mass of 0.300 keV. This value of the pion 
momentum has been choosen because, at the JINR Phasotron, the DUBTO 
experiment [13] studying pion interactions with 4 He nuclei at 200 MeV/c, 
collected several 7r— events. In this case the allowed values of the <fi an- 
gle are less than 11 deg. Setting the value of the angle <fi equal to 5 deg, 
we applied the procedure described in paragraph 3.2. The uncertainties of 
the momenta were varied in the range 10 2 -^10~ 9 MeV/c and of the angle <fi 
in the range 1-^10~ 5 deg. Figure 5 shows the same m v , m* and m** upper 
limits as in Figure 4, calculated using the m 2 ,, m* and m** distributions as 
functions of the momentum uncertainties in the range 10 2 -i-10~ 9 MeV/c, and 
with an error for the angle <f> of 10~ 4 deg. As one can see, at 200 MeV/c the 
m v upper limit values follow the line (a) of Figure 4, which was drawn for 
pions decaying at lower energies. The behaviours of the m* and m** upper 
limits appear to be nearly parallel to those of Figure 4. 

Considering the experimental uncertainties actually attainable in labora- 
tories and taking into account the present values of pion and muon masses 
and their uncertainties, it appears impossible to perform experiments, based 
on the measurement of all the kinematic parameters of the 7r^/iz/^ decay, 
from which an upper limit of the muon neutrino mass less than about 1 keV 
can be deduced. 

In the frame of the results of the present analysis, using the distributions 
of Figure 3 relative to the measurement of the event displayed in Figure 2 
and with the uncertainties of the present experiment, we obtained the upper 
limits of m„ and m*, at given confidence levels, reported in Table 2 using 
the classical statistical approach. The values have been deduced with formu- 
lae (7); in particular, m* = (-0.186+1.282*0.211) MeV = 0.0845 MeV, at 
90% confidence level. It must be noted that, with the upper limit of m** = 
59.96 MeV, one obtains 

Vm*m** = V0.0845 • 59.96 MeV = 2.2 MeV. (11) 

The above values of m*, m** and m„ upper limits are reported in Figure 4. 

4 Conclusions 

This paper describes the determination of a confidence upper limit for the 
muon neutrino mass by measuring the momenta of the decaying pion and of 
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Figure 5: Behaviour of the upper limits of m v (open stars), m* (full stars) 
and m** (open rhombs) calculated from the ml, m* and m** distributions at 
90% confidence level, as a function of the experimental uncertainties in the 
momenta, and for A<p = 10~ 4 deg, for pions of 200 MeV/c. The lines (a) and 
(b) are those drawn to guide the eye in Figure 4. The line (c) represents the 
upper limit value of 60 MeV. 
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m u (MeV) 


m* (keV) 


p(%) 


1.4 


33 


80 


2.2 


85 


90 


3.1 


161 


95 



Table 2: Upper limits of the muon-neutrino mass (deduced from m%) and of 
m* at a p% confidence level, for the event of Figure 2 measured and analysed 
by the method described in the text. 

the produced muon in a unique (nfxe) nuclear event recorded in the PS 179 
experiment at LEAR of CERN. A sizeable improvement of previous results 
obtained using low density visualizing detectors has been achieved. 
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